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Aims. To find out whether toroidal field can stably exist in galaxies the current-driven instability of toroidal magnetic fields is consid- 
ered under the influence of an axial magnetic field component and under the influence of both rigid and differential rotation. 
Methods. The MHD equations are solved in a simplified model with cylindric geometry. We assume the axial field as uniform and 
the fluid as incompressible. 

Results. The stability of a toroidal magnetic field is strongly influenced by uniform axial magnetic fields. If both field components 
are of the same order of magnitude then the instability is slightly supported and modes with m > 1 dominate. If the axial field even 
dominates the most unstable modes have again in > 1 but the field is strongly stabilized. All modes are suppressed by a fast rigid 
rotation where the m = 1 mode maximally resists. Just this mode becomes best re-animated for CI > Q.a (£2a the Alfven frequency) if 
the rotation has a negative shear. - Strong indication has been found for a stabilization of the nonaxisymmetric modes for fluids with 
small magnetic Prandtl number if they are unstable for Pm = 1. 

Conclusions. For rotating fluids the higher modes with m > 1 do not play an important role in the linear theory. In the light of our 
results galactic fields should be marginally unstable against perturbations with m < 1. The corresponding growth rates are of the order 
of the rotation period of the inner part of the galaxy. 
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1. Introduction 

It seems to be an open question whether or not toroidal mag- 
netic fields are stabilized under the presence of poloidal field 
components. Lundquist (1951) formulated that azimuthal fields 
with energy exceeding the energy of the axial field component 
become unstable. With other words, he found that uniform ax- 
ial fields act stabilizing what - if true - would form a highly 
interesting finding also for MHD experiments in the laboratory. 

By use of a cylindric magnetic geometry Roberts (1956) 
opened a new discussion and found for all ratios of the az- 
imuthal field component and the axial field component instabil- 
ity against perturbations with high azimuthal mode numbers m. 
Tayler (1980) re-discussed the adiabatic stability of stars with 
mixed poloidal and toroidal fields. For poloidal and toroidal field 
components of the same order he suggested stability of the sys- 
tem but the final answer to the question remained open. In his 
detailed paper about magnetic instabilities Acheson (1978) only 
considered the stability of purely toroidal fields. Extending this 
work one can ask for the stability of a simplified magnetic model 
where axisymmetric and stationary toroidal field are considered 
under the influence of a homogeneous axial field which itself 
is stable by definition. For an ideal medium Bonanno & Urpin 
(2010) considered the stability of such constellation without any 
rotation with respect to applications in jet theory. They exclude 
stability for fields with pitch |B Z |/|B^| of order unity. Particular 
attention is given to the instability of nonaxisymmetric modes 
with azimuthal mode numbers m > 1. If the axial field domi- 
nates, the instability persists for rather high mode numbers (of 
order 100). The latter has been found by use of an almost iden- 



tical model by Tayler (1960) including the differences of the so- 
lution for different helicity of the background field. 

A more complicated model has been investigated by 
Braithwaite (2009) where also the poloidal field component (ax- 
isymmetric as also the toroidal one) can be unstable alone but the 
author finds stability of the combination of poloidal and toroidal 
field if both components are of the same order. 

We shall show for dissipative fluids that compared to the 
case of purely toroidal fields the configuration with \B Z \ ^ \B^,\ 
is (slightly) more unstable while increasing \B Z \ more and more 
stabilizes the toroidal fields. With rotation the situation changes. 
If the rotation rate exceeds the Alfven frequency of the toroidal 
field the instability of any m is suppressed but the mode with 
m = 1 persists at longest. If, as it is almost always the case, 
the rotation is not rigid then only the modes with small m sur- 
vive and the axisymmetric standard MRI starts to dominate for 
sufficiently fast rotation. We find that the competition between 
the modes m = and m = 1 should be observable with galax- 
ies. Note that only a few examples of nonaxisymmetric mag- 
netic field patterns have been found by the observers (Beck et al. 
1996). 

Probably, the hydromagnetic jets are suitable subjects for 
the application of magnetic field instabilities but for the low 
Reynolds numbers and Hartmann numbers which we shall deal 
with in the present paper the galaxies containing a remarkable 
interstellar turbulence are forming better objects. Galaxies pos- 
sess quadrupolar-type magnetic fields with toroidal and poloidal 
components of the same order of magnitude (~ 10~ 5 Gauss) and 
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with the phase relation BrB^ < 0. They rotate with the charac- 
teristic rotation law 



Q.{R) 



const. 
R 



(1) 



with RQ. ^ 200 km/s. Their characteristic density is of order 
10~ 24 g/cm 3 and the magnetic diffusivity is about 10 26 cm 2 /s due 
to the action of interstellar turbulence. The resulting magnetic 
Reynolds number Rm ^ URo/rj is thus of order 1000 while the 
Lundquist number of the toroidal field S = BRq/ yj^pr] reaches 
values of 200. For galaxies both characteristic numbers are thus 
of the same order or more precisely 



0^5 0. 



(2) 



expressed with the rotation rate and the Alfven frequency Qa 
(see Eq. ( fTol l. below). The question is whether such a mag- 
netic constellation is stable against nonaxisymmetric distur- 
bances with the azimuthal mode number m. 

The stability of a toroidal field strongly depends on its radial 
profile. So the current-free profile oc l/R is stable against 
disturbances with the azimuthal mode numbers m = 0,1,.... 
On the other hand, the profile B$ oc R is stable only against 
m = but it is unstable against disturbances with m > (Tayler 
1957; Velikhov 1959). The latter profile will mainly be used in 
the present paper by its simplicity as it is due to a homogeneous 
axial electric current. A sufficiently strong toroidal field which is 
nearly uniform in radial direction is also unstable against distur- 
bances with the mode number m > 0. We shall model a galactic 
magnetic field with respect to the equator by means of a Taylor- 
Couette flow with the rotation law (Q]) periodic in the axial di- 
rection z. The axial electric current which produces the toroidal 
magnetic field component is assumed as homogeneous in z. It 
is obvious that such a simple cylindric model cannot describe 
the field geometry in a global rotating disk but only in one of its 
hemispheres. 

For simplicity the cylinders which confine the conducting 
fluid are highly conducting, and no-slip boundary conditions are 
used at the cylinder walls. The magnetic background field (as- 
sumed as stationary) also possesses a uniform axial field com- 
ponent so that the resulting field pattern forms a spiral. It is the 
stability of such a spiral with fixed current helicity which is con- 
sidered in the present paper. With respect to galactic applications 
this is an over-simplification as for dynamo-generated magnetic 
fields of either parity (with respect to the equator) the current 
helicity always behaves antisymmetric. 

It appears to be reasonable not to limit the azimuthal mode 
number to \m\ 5 1 so that also higher values can be considered 
(see Arlt et al. 2007; Bonanno & Urpin 2010). 



The resulting linear equations are 

^ + (U ■ V)« + (« • V)17 = --Vp + vAu + 
at p 

H — — curl b x B + — !— curl B xb, 

y"0P y"0P 



db 
dt ' 
and 



curl(n x B) + curl(t7 x b) + T]Ab, 



div u = div b — 0, 



(5) 
(6) 

(7) 
(8) 



where u is the perturbed velocity, b the perturbed magnetic field, 
p the pressure perturbation and v and rj are the kinematic viscos- 
ity and the magnetic diffusivity. The magnetic Prandtl number is 
defined by 



Pm = -. 



The stationary background solution is 



b B 

Bq, = agR + — , 
K 



(9) 



(10) 



where an, bn, cib and bg are constants defined by 



an = £2 



m-w 



a B = 



1-fj 2 

Bin flil^B ~ fj) 



Rin 1 



bn = Q,r,R;. 



b B = B m Ri 
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(11) 



(12) 



Here Ri n and R out are the radii of the inner and outer cylinders, 
fii n and Q out are their rotation rates, and B m and B out the az- 
imuthal magnetic fields at the inner and outer cylinders, resp. 

The outer value B out is normalized with the uniform vertical 
field, i.e. 



(13) 



As usual, the toroidal field amplitude is measured by the 
Hartmann number 



Ha = 



VWV77' 



(14) 



2. The equations 

We are interested in the linear stability of the background field 
B = (0, B^R), B ), with B = const, and the flow U = 
(0, RQ.(R), 0). The perturbed system is described by 



u R , m , u z , p, b R , b^, b. 



Developing the disturbances into normal modes, the solutions of 
the linearized MHD equations are considered in the form 



/ = /CR)exp(i(fe + m<p + art)), 



(4) 



where / is any of the velocity, pressure, or magnetic field distur- 
bances. 



Here Rq = V^in(^out - ^in) is used as the unit of length, r]/Ro 
as the unit of velocity and B m as the unit of the azimuthal fields. 
Frequencies, including the rotation £2, are normalized with the 
inner rotation rate Qj n . The ordinary Reynolds number Re and 
the magnetic Reynolds number Rm are defined as 



(3) Re 



Qd2 
m 



Rm 



(15) 



and the Lundquist number S is defined by S = Ha ■ VPm so 
that it can also be understood as the magnetic Reynolds number 
formed with the Alfven frequency, 



B„ 



(16) 



G. Riidiger et al.: Current-driven instability of rotating helical magnetic fields 



3 



instead of the rate of the global rotation. Note that the normal- 
ization concerns the maximal value of the toroidal field. 

The boundary conditions associated with the perturbation 
equations are no-slip for u, i.e. ur — = u z — 0, and perfectly 
conducting for b, i.e. db^/dR + b$/R — bR — 0. These boundary 
conditions hold for both the inner and the outer cylinder. 

All our calculations refer to a container with R out = 2R in , 
i.e. fj - 0.5. For this choice a field which is current-free in the 
fluid is described by yUg = 0.5. A homogeneous axial electric 
current between the cylinders requires fig = 2 which is the pre- 
ferred value in this paper. The axial magnetic field component 
is assumed as uniform so that the resulting current helicity of 
the magnetic field is also homogeneous. This can only be true 
within one hemisphere of the celestial body (here galaxy). Our 
model does not define an equator. We shall consider that hemi- 
sphere where the current helicity B ■ J of the background field 
is negative, i.e. B < 0, which is an arbitrary choice. As shown by 
Riidiger et al. (2010) the resulting instability forms left spirals if 
the rotation is slow. The kinetic helicity (u ■ curl w) of the pertur- 
bations (averaged over the azimuth) proves to be positive for this 
field. The kinetic helicity does not change its sign if the rotation 
is faster but then the magnetic pattern forms right spirals. 



3. No rotation 

We start to consider a nonrotating container, i.e. Re = 0. In 
this case for given geometry and given vector (O^^R), Bo) of 
the magnetic field the critical Hartmann number does not de- 
pend on the magnetic Prandtl number Pm. The azimuthal drift 
of the nonaxisymmetric instability pattern vanishes (see Riidiger 
& Schultz 2010). Without rotation the instability patterns do not 
drift in azimuthal direction. For rapid rotation the drift rate al- 
ways grows with Qj n rather than with B^ (see below). 
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Fig. 1. Neutral stability for negative/? and for/z^ = 2. The curves 
are marked with their mode number m. The minimum Hartmann 
number of the toroidal field with Ha = 57.6 exists for \B\ - 8 
(thin horizontal line). The plot valids for all magnetic Prandtl 
numbers. 

The most natural stationary magnetic profile is B^, cc R which 
is the result of a homogeneous electric current flowing through 
the whole domain with R < R out (see Roberts 1956). We know 
that for this case and for very large \B\ the critical Hartmann num- 
ber has the value of 70.6 for m — 1 (see Riidiger et al. 2007). 
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Fig. 2. The wave numbers normalized after ([TBI for various B. 
Note that the physical wave numbers run with the value of \B\. 
fig = 2, Pm = 1. 

This value of the critical Hartmann number is slightly re- 
duced if a small and uniform axial component of the magnetic 
field is added to the system. Hence, a uniform axial field supports 
the pinch-type instability of the toroidal field. This effect, how- 
ever, is rather weak: the critical Hartmann number sinks from 
about 70 to about 60 (see the horizontal line in Fig.Q]). For m > 1 
the destabilization of the toroidal field by axial fields is much 
stronger so that for \B\ of order unity all the modes with different 
m have more or less the same critical Ha. We thus do not find 
a stabilizing effect to toroidal fields by axial fields components 
compared to fields of purely toroidal fields. 

For B - -8 we find Ha = 57.5 as the absolute minimum 
of the stability curve for m = 1. For stronger axial fields the 
critical Hartmann number basically grows reaching values of 
about 1000 for B -0.1. For strong axial fields the modes 
m > 1 possess lower critical Hartmann numbers than the mode 
with m = 1. The differences of the curves with various m are 
much smaller than those for weak Bo but the Fourier compo- 
nent with m — 4 possesses the lowest critical Hartmann num- 
ber for B - -0. 1 . Nevertheless, for dominating axial field the 
toroidal field is strongly stabilized - in particular the mode with 
m — 1 . The modes with m > 1 possess somewhat slower critical 
Hartmann numbers but also these modes are basically stabilized 
(see Fig.[T|). 

To summarize: The pinch-type instability of toroidal fields 
under the presence of a uniform axial magnetic field without ro- 
tation is strongly suppressed by strong axial fields. The maxi- 
mal stabilization happens for m = 1, hence the most unstable 
modes have azimuthal mode numbers m > 1. If B^ and B z are 
of the same order then the field is (slightly) more unstable than 
for B z = 0. We find that with strong enough axial current-free 
magnetic fields rather strong toroidal fields can be stored in the 
container. 

All results are invariant against the simultaneous transforma- 
tion m — > — m and B — > -B. 

The wave numbers k of the unstable modes reflect the insta- 
bility pattern. The shape of the cells is described by the relation 



Sz 



Rnut — ^ii 



n 
k' 



(17) 



The critical wave number for purely toroidal fields with hb — 2 
is 2.8 so that after (TTTb the cells are almost spherical. The wave 
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numbers with axial field component are shown in Fig. [2] As ex- 
pected they linearly grow for growing |m| and for growing 
For dominating axial field the cells become longer and longer 
what is not unexpected. For B < - 1 the simple relation 



— -0.5, 
mB 



(18) 



results so that the aspect ratio of the cells is more and more given 
by the pitch of the field, i.e. 



Sz 
i — 
6R 



B : 
B,i, 



(19) 



After Fig.|2]this relation is well-established for \B Z \ > \B$ 
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Fig. 3. No rotation: The growth rate in units of the diffusion fre- 
quency runs linearly with the Hartmann number of the toroidal 
field. Top: B = -20. The mode with m — 2 has the largest growth 
rate. Bottom: B — -2 (see the vertical lines in Figs. Q] and [2]). The 
curves are marked with their mode number m. Generally, for fi- 
nite B maximum growth rates occur for m > 1. fis = 2, Pm = 1. 

The growth rates y = -3(w) must be given in units of the 
diffusion frequency J]/R^. We shall find that at least the modes 
with the higher m are strongly dependent on the Reynolds num- 
ber of rotation. Without rotation the growth rates for given B and 
Pm are plotted in Fig. [3] The used B close to the minimum where 
\B<f\ - \B Z \ will be the preferred value for many of the examples 
presented in this Paper. 



Both plots in Fig. [3] demonstrate the finding that the growth 
rates run with the magnetic Alfven frequency Qa- For stronger 
fields strong differences for the growth rates of various m appear. 
For dominating azimuthal field (\B\ = 20, top) this is a weak 
effect but for and B z of the same order (\B\ = 2, bottom) it 
is strong. Of course, there are maxima; but for higher Hartmann 
number the highest growth rates belong to higher m. 

The dependence of the growth rates on the magnetic Prandtl 
number is a complex problem. The majority of the numerical 
simulations concerns to Pm = 1. In the Sect. 5.1 below we shall 
show that indeed this choice forms a special case. For resting 
cylinders the nonaxisymmetric mode with m — 1 grows fastest 
for Pm = 1 if it is normalized with the geometrical average 
rf = yjvrj of both diffusivities. For given product of v and r\ 
the mode for almost equal diffusivities is most unstable while it 
becomes even stabilized if the two viscosities are too different. 
The consequences for this finding may be dramatic if applied 
to numerical simulations. A field may be unstable for Pm = 1 
which proves to be stable for more realistic very small or very 
large Pm. 



4. Rigid rotation 

It is known that rigid rotation stabilizes the magnetic perturba- 
tions. This effect can easily be realized with our model. For the 
standard model with fig - 2 and Pm = 1 the growth rates have 
been calculated together with the drift rates for a supercritical 
value of Ha. 

We start with a very small pitch, i.e. with nearly toroidal 
fields (B = -20). Figure [4] gives the results for the growth rates 
normalized with the diffusion frequency. Both the given modes 
for m = 1 and m — 2 are strongly suppressed by the basic ro- 
tation. We also find, however, that the mode with m = 1 better 
survives the rotational suppression than the modes with higher 
m. 
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Fig. 4. Rigid rotation and nearly toroidal fields (B = -20): The 
growth rates divided by the diffusion frequency for supercritical 
magnetic field with Ha = 200. The magnetic instability is sta- 
bilized for Re > Ha. The kink-type mode with m = 1 survives 
rotation better than the higher modes, fig — 2, B — -20, Pm = 1. 



After Fig. Q] the most interesting situation should exist for 
magnetic fields with a pitch angle of order unity. The eigenvalues 
for the field with B — -2 have thus been calculated. Figure 
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gives the main results. The critical Ha for /3 — -2 after Fig. Q] 
is ~ 70. One finds positive growth rates for slow rotation and 
stability for fast rotation. The instability cannot exist for Re > 
Ha or in other words, for 

Q > Q A . (20) 

Again the mode with m = 1 withstands at best the rotational 
suppression. It is also true that the modes with the highest m are 
suppressed already by lower Reynolds numbers. Obviously, the 
dominance of the modes with m > 1 disappears by rigid rotation. 

For fast rotation the drift !R(w)/mQ of all modes approaches 
the value -1 so that after the relation 

df m 

an observer in the laboratory system finds the magnetic pattern 
as almost corotating. If the stellar rotation can only be observed 
via their magnetic variation then the rotation of such an object is 
well approached by the rotation of the magnetic pattern. The ob- 
served (magnetic) rotation is slightly slower than the real plasma 
rotation. However, there is a jump in the curves: for slow rota- 
tion the modes with m = 1 and m — 2 rotate much faster than the 
container. We have found such a jump already for rotating stars 
with unstable toroidal fields (Riidiger & Kitchatinov 2010). 

5. Differential rotation 

There is a very new situation if the outer cylinder rotates slower 
than the inner one. The simplified rotation law of our model may 
be the galactic one, ([]]), so that for fj = 0.5 the rotation ratio 
is fin = 0.5 which also can be considered as the normalized 
rotation of the outer surface of the container. 

The growth rates for Ha = 160 and the mentioned differen- 
tial rotation are given in Fig. [6] (top). The plot is identical with 
the plot for rigid rotation (Fig. [5] top) if the rotation is slow, i.e. 
for Q. <k Q.A- All modes are rotationally stabilized. For Q > 
however, the magnetic instability is re-animated but at most for 
the lower modes. Finally the m = 1 mode becomes dominant; its 
growth rate (in diffusion units) becomes higher and higher finally 
running with the rotation frequency. This new type of magnetic 
instability even exists for current-free toroidal magnetic fields 
so that we have named it the azimuthal magnetorotational insta- 
bility (AMRI). It is basically nonaxisymmetric with low m and 
results from the interaction of differential rotation and toroidal 
fields (Riidiger & Schultz 2010). The growth rate y runs with Q 
rather than Qa when Q. > Qa (Fig.[6]l. 

The higher modes dominate only for small Reynolds num- 
bers. They do not contribute to the instability for high Reynolds 
numbers as they are damped by fast differential rotation. As we 
have shown in Sect. 3 the modes with m > 1 are also damped 
for weak and for strong extra axial magnetic field components. 
Where they are most unstable (for/? of order unity) any rotation 
does suppress them. Their domain of dominance should thus be 
small in astrophysical applications. 

5.1. Standard MRI 

Also axial fields can be unstable under the presence of differ- 
ential rotation. The leading mode of this standard MRI is ax- 
isymmetric (see Kitchatinov & Riidiger 2009). We have given in 
Fig- also the growth rate of the mode m = 0. It possesses the 
largest growth rate y if the rotation rate Q is high enough, in this 
case the growth rate also runs with Q. We find that for D. » 
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Fig. 5. Rigid rotation (pn = 1): Top: Growth rates in unit of 
diffusion frequency. The magnetic instability is stabilized for 
Q. > Oa- Bottom: the azimuthal drift after (l2"TT i of the modes. 
p B = 2,/3 = -2, Ha= 160, Pm = 1. 



the most interesting case of \B$\ — \B Z \ leads to a dominance of 
the standard MRI. Note the existence of an intersection between 
the growth rates of m = and m = I. Left from this point the 
nonaxisymmetric mode dominates the axisymmetric one while 
right of this point it is opposite. After Eq. (|2]i galaxies do exist 
very close to that point. One should thus be aware that the sta- 
bility of galactic fields should be rather delicate. We have shown 
that the negative shear of the rotation law strongly destabilizes 
the toroidal field. It is thus not clear, however, whether the most 
unstable mode is axisymmetric or not. 

5.2. The Pm-dependence of the growth rates 

The dependence of the growth rates on the magnetic Prandtl 
number Pm is not trivial. We have shown that for resting contain- 
ers the characteristic Hartmann numbers for marginal instability 
do not depend on Pm (Riidiger & Schultz 2010). This is not true 
for the growth rates - and this the more the faster the rotation 
is. Additionally, it is not obvious how to normalize the growth 
rates and the Reynolds numbers if both the diffusion times dif- 
fer. We have also shown that the use of frequencies normalized 
with a geometrically averaged diffusion 77* which is symmetri- 
cally formed with v and r\ with rf - ^Jvrj = r\ VPm seems to be 
most appropriate. 
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Fig. 6. Differential rotation (fin = 0.5): Growth rates in unit 
of diffusion frequency. The kink-type instability (m = 1) is re- 
animated by fast rotation (D. > Q.a)- Hb — 2, /3 — -2, Ha = 160, 
Pm=l. 

In Fig. Qthe growth rates without and with (differential) ro- 
tation are given for a fixed Hartmann number. Both the growth 
rates of the mode m = 1 and the global rotation rate are normal- 
ized with rf, hence it is 



One finds Pm = 1 always leading to maximum growth rates 
for slow and fast rotation. Either small or large magnetic Prandtl 
numbers lead to slower growth of the instability than for Pm = 1 . 
This effect is so strong that the considered field pattern can even 
be stabilized if the magnetic Prandtl number is too small or too 
high. This is indeed the case in Fig. [7] for Pm < 0.01. An in- 
stability found with numerical simulations for Pm = 1 does not 
automatically exist for much smaller or much larger Pm. If for a 
given value of rj* the magnetic field is unstable for Pm 1 this 
must not be true if the numerical values of v and 77 are too differ- 
ent. This is an important restriction of the validity of numerical 
simulations of magnetic stability/instability which are operating 
with v = tj. The stability/instability of magnetic fields strongly 
depends on the magnetic Prandtl number of the fluid. For resting 
or rotating media the fields are most unstable for Pm = 1 . 

Note that for given Hartmann number (here Ha = 160) one 
finds in Fig. [7] two regimes for the rotational influence on the 
growth rates. There is almost no influence of small Rm* on the 
growth rate a>* . Figure [7] (top) shows the maximum growth rate 
(for Pm = 1) as of order 10 leading to a minimum growth time 
of 0.1 diffusion times. For galaxies with Rq =T0 kpc and with 
tj m 10 26 cm 2 /s the diffusion time is then 3 Gyr. One finds, how- 
ever, fast global rotation accelerating the instability. From Fig. 
[7] (bottom) the physical growth rate results as Q/5 so that the 
growth time is reduced by the rotation to about one rotation time. 
The rotation of the inner part of the galaxy is here concerned, 
with rotation times of about 50 Myr. Hence, the current-driven 
magnetic instability is a rather fast process. 

6. Almost homogeneous toroidal field 

To check the consistency of our model in more detail the Figs. 
Q]and[6]are modified for almost uniform toroidal fields, i.e. with 




250 500 750 1000 
Re* 

Fig. 7. The growth rates for the modes m = 1 normalized with if 
(see text) for fixed Hartmann number and for various Pm. Top: 
no rotation, Bottom: quasi-galactic differential rotation (//n = 
0.5). The curves are marked with their magnetic Prandtl number. 
Ha= 160,/?= -2, (i B = 2. 



Hb — 1. We find very similar results but with slight numerical 
differences. Figure [8]reveals the stabilizing action of axial fields 
to the pinch-type instability of the toroidal field to be much more 
effective for the case of almost homogeneous B^. The critical 
Hartmann number for instability grows by orders of magnitudes 
if the axial field grows only by a factor of five. 

Also the complex influence of differential rotation on the in- 
stability of those fields with same order of toroidal and poloidal 
magnitude shown by Fig. |6]exists for the case of almost homo- 
geneous toroidal fields (Fig.|9]l. Slow rotation acts stabilizing but 
the modes with higher m dominate for a while. For fast rotation 
the modes with low m are re-animated as they form the new in- 
stability. For slow rotation the modes with higher m exhibit the 
maximum growth rates but for fast and differential rotation the 
mode with m = 1 grows fastest. 

Again the axisymmetric mode with m = is concerned 
which starts to dominate beyond the given crossing point. The 
coordinates of the intersection are nearly the same as in Fig. [6] 

Obviously, the above findings about the influence of differen- 
tial rotation on the stability of magnetic fields with spiral struc- 
ture do not basically depend on the radial profile of the toroidal 
field. It makes thus sense to call - as we did - this effect as the 



G. Riidiger et al.: Current-driven instability of rotating helical magnetic fields 



7 



10000 

Ld 
CD 

I 1000 



00 



a: 
< 
x 



\\ 1 


Re = 








1 



-0.1 



-1 

|8 



-10 



Fig. 8. The same as in Fig. Q]but for fig = 1- The absolute 
minimum of the Hartmann number of Ha = 86.8 is located at 
B — 7.25 for m — 1. Note the extreme stabilization of for 
increasing B z . 

magnetorotational instability which originally only concerned 
current-free toroidal fields (fig = 0.5, in our notation). 
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Fig. 9. The same as in Fig.|6]but for/i# = 1. The dashed line is 
for m = which dominates for faster rotation. With the values 
characteristic for galaxies (right axis) one finds that the axisym- 
metric standard MRI should dominate. 



7. Summary 

In a cylindric geometry the pinch-type instability of axisymmet- 
ric and also magnetic spirals with finite current helicity B ■ J is 
considered under the influence of rotation. The field is formed by 
an unstable toroidal field and a uniform axial field which is sta- 
ble by definition. The pitch of the spiral is given by the inverse 
of B = B t j,IB z which is negative for the considered lefthanded 
spirals. The larger the pitch of the background field the higher 
is the azimuthal Fourier m of the mode with the largest growth 
rate. 

The excitation of modes with low m is here of particular in- 
terest. For small \B\ and without rotation typically a mode with 
m > 1 is excited with the largest growth rate. As stressed by 



Bonanno & Urpin (2010) this phenomenon could have conse- 
quences for the jet theory. As we have shown, however, a grow- 
ing axial field stabilizes the toroidal field more and more. The 
critical Hartmann number grows by orders of magnitudes if 
\B\ reduces from order unity to order 0.1 (see Figs. [TJ and [SJ. 
Helical background fields with large axial field component are 
thus much more stable than purely toroidal fields without finite 

Also a global rotation stabilizes the pinch-type instability. 
Figure|5]shows for a magnetic field with almost equal field com- 
ponents (B = -2) how the rotation quickly stabilizes the modes 
with m > 1 while the kink-instability (m - 1) remains unstable 
for a little faster rotation. The growth rates of the modes are con- 
tinuously reduced by growing Reynolds numbers. Generally, the 
helical background fields are stable against all nonaxisymmetric 
perturbations if Q. » Qa- 

A very new situation results for nonrigid rotation. Figure [6] 
clearly demonstrates with a rotation law known from galaxies 
that for Q > Qa the growth rates after a characteristic minimum 
at Q ^ Qa again reach positive and large values. A rotation 
law with positive shear will always stabilize the nonaxisymmet- 
ric instability. It has also been shown that under the presence 
of differential rotation with negative shear the toroidal field can 
become unstable even if there is no electric current in the con- 
tainer (see Riidiger & Schultz 2010). While for slow rotation the 
modes with higher m are most unstable it is for fast rotation the 
mode with m = I. 

If the field possesses an axial component then under the in- 
fluence of differential rotation with negative shear the standard 
MRI appears in form of an growing axisymmetric (m = 0) roll. 
The lines of marginal instability for m = and m — 1 are cross- 
ing so that for fast enough rotation the axisymmetric perturba- 
tion dominates. In any case we find that a spiralic magnetic field 
under the influence of differential rotation with negative shear 
appears to be extremely unstable (see Fig. [9] for a field nearly 
uniform in radius). As our model roughly reflects the magnetic 
geometry in galaxies with SN-driven interstellar turbulence one 
should expect their dynamo-generated magnetic fields as rather 
unstable. Nonaxisymmetric global field configuration might not 
be the exception. 

Our model also allows the variation of the magnetic Prandtl 
number. Figure [7] shows a rather clear situation. Its parameters 
only depend on the product of nu and 77, i.e. they are invariant 
against an exchange of v and 77. We find for fixed value of v • 77 
strong differences of the lines for v = 77 and v + 77. The me- 
dia with v = 77 are more unstable than the media with v + 77. 
Moreover, if the magnetic field is unstable for Pm = 1 it can 
even be stable for Pm » 1 or Pm <K 1. The result requires 
care with the interpretation of numerical instability calculations 
if the considered medium has a magnetic Prandtl number much 
smaller than unity. A magnetic field configuration which for a 
given Hartmann number and Pm = 1 results as unstable can be 
stable for much smaller or much larger Pm. 
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